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Abstract. Single- particle transport in disordered potentials is investigated on scales 
below the localization length. The dynamics on those scales is concretely analyzed 
for the 3-dimensional Anderson model with Gaussian on-site disorder. This analysis 
particularly includes the dependence of characteristic transport quantities on the 
amount of disorder and the energy interval, e.g., the mean free path which separates 
ballistic and diffusive transport regimes. For these regimes mean velocities, respectively 
diffusion constants are quantitatively given. By the use of the Boltzmann equation 
in the limit of weak disorder we reveal the known energy-dependencies of transport 
quantities. By an application of the time-convolutionless (TCL) projection operator 
technique in the limit of strong disorder we hnd evidence for much less pronounced 
energy dependencies. All our results are partially confirmed by the numerically exact 
solution of the time-dependent Schrodinger equation or by approximative numerical 
integrators. A comparison with other findings in the literature is additionally provided. 
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1. Introduction 

Any solid contains disorder: Either there are impurities, vacancies, and dislocations in 
an otherwise ideal crystal lattice. Or there is no lattice structure at all. An abstract 
quantum system which is commonly used as a paradigm for transport in real disordered 
solids is the Anderson model pQ . In its probably simplest form without particle-particle 
interactions (electron-electron, electron-phonon, etc.) the Hamiltonian may be written 
as 



where a? and a? denote the usual annihilation, respectively creation operators; f labels 
the sites of a d- dimensional (cubic) lattice; NN indicates a sum over nearest neighbors 
r*i and r* 2 ; and represent independent random numbers, e.g., according to a Gaussian 
distribution with mean (e?) = and variance (e^ e? 2 ) = Sf lt p 2 a 2 . Even though such a 
distribution is considered throughout this work, the random numbers can be realized 
according to a Lorentzian, box, or binary distribution as well [2j [31 Hj. In all cases 
disorder is implemented in terms of a random on-site potential. (Random hopping 
coefficients are sometimes taken into account, too.) 

In the presence of such a disorder, cr^O, the eigenstates of the Hamiltonian are no longer 
given by Bloch functions: Instead the eigenstates are not necessarily extended over the 
whole lattice and can become localized in configuration space, i.e., the envelope of a 
wavefunction decays exponentially on a finite localization length (2j |5]. The finiteness 
of the localization length is one manifestation and, say, definition of the localization 
phenomenom. (There certainly are other mathematical definitions of localization, e.g., 
the finiteness of the inverse participation number, the independence of eigenvalues from 
boundary conditions, etc. [2]) This phenomenon and particularly its impact on transport 
have intensively been studied for the Anderson model [Il[2l|3jlIJ|5j[6j[7JIEllHl ^ ' 
For the lower dimensional cases, d = 1 and d = 2, all eigenstates of the Hamiltonian 
feature finite localization lengths for arbitrary (non-zero) values of o (except for 



e ? d\d ? + ^2 
NN 



(1) 
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situations with short-range correlated disorder, e.g., as realized in the random dimer 
model [IIIEG])- Therefore in the thermodynamic limit, i.e., with respect to the infinite 
length scale an insulator is to be expected. Of particular interest is the 3-dimensional 
case, as considered in the work at hand. Here, a mobility edge, i.e., a certain cross-over 
energy separates the spatially localized from the spatially extended wavefunctions in 
energy space [21 El [5]. When the amount of disorder is increased, the mobility edge goes 
above the Fermi level and a metal-to-insulator transition is induced at zero temperature, 
still with respect to the infinite length scale. When a is further increased above the 
critical disorder {p c ^ 6 for a Gaussian distribution [21 S]), all eigenstates become 
localized and an insulator is to be expected for each temperature (without particle- 
particle interactions). The Anderson metal-to-insulator transition is widely believed to 
be continuous without a minimum conductivity, e.g., as supported by the one-parameter 
scaling theory of localization [21 El [7] . 

Our work, other than most of the pertinent literature, focuses on the dynamics on 
scales below the localization length. We particularly intend to analyze the dynamics on 
those scales comprehensively as a function of energy and disorder. Qualitatively, our 
analysis allows to identify two regimes of length scales which are purely ballistic and 
strictly diffusive (rather than super diffusive, subdiffusive, or anything else). It generally 
is a challenge to theoretically confirm reliably the "presence of diffusion" in strongly 
disordered and/or interacting quantum systems. Quantitatively, our analysis enables 
the evaluation of mean velocities, respectively diffusion coefficients for a wide range of 
disorders between zero and the vicinity of oc- Such a detailed knowledge about diffusion 
constants appears to be important, especially since dc-conductivities are directly related 
by the Einstein relation, at least for a < o~c- In the limit of strong disorder diffusion 
coefficients have been suggested in the literature by the numerical study of Green's 
functions for very few disorders and a single energy at the spectral middle solely [H [9] . 
The dependence of diffusion constants on energy is usually discussed in the limit of weak 
disorder only. However, also in that limit, we demonstrate that energy dependencies are 
much richer than common approximations for a free electron gas [2]. 
The work at hand is structured as follows: First of all we provide a qualitative picture of 
the dynamics on scales below the localization length in section [21 Then this qualitative 
picture is subsequently developed and quantitatively confirmed in the whole section |3l 
The limit of weak disorder is firstly analyzed in section 13.11 by the use of the Boltzmann 
equation [131 ttH H3 EE]. The limit of strong disorder is afterwards investigated by 
an application of a method which is mainly based on the time-convolutionless (TCL) 
projection operator technique [TFJ [TH1 HH [2Ql EI]. In section 13.21 the method as such 
is introduced and its predictions on the dynamics are presented. The validity range 
of these predictions is analytically discussed in section 13.31 and numerically verified in 
section 13.41 We finally close with a summary and conclusion in section HI 
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2. Qualitative picture of the dynamics 

In the present section we intend to provide firstly a qualitative picture of the dynamics 
on scales below the localization length. This qualitative picture essentially summarizes 
the findings of the methods which are introduced in detail and applied concretely in 
the following sections. In particular we emphasize the main conclusions of the work at 
hand and discuss these conclusions in the context of known results in the literature. In 
this way we also give a comprehensive summary for the readers which are not primarily 
interested in the methodic details. Apart from that the summary certainly makes the 
line of thoughts in the subsequent sections more plainly. 

The above mentioned qualitative picture of the dynamics is illustrated in figure HJ In 



(a) weak disorder (b) strong disorder 




energy E 



Figure 1. Sketch for the dependence of transport on the wave number q and the 
energy E for (a) weak disorder and (b) strong disorder. Both top panels indicate 
the rough position of localized (loc, gray), diffusive (dif, white) and ballistic (bal., 
black) regimes. Particularly, a possible "corridor" of wave numbers is indicated (dashed 
lines), where transport is diffusive for almost all energies. Both bottom panels display 
suggestions for qualitatively different energy dependencies of the diffusion constant. 

this figure a sketch for the dependence of transport on the length scale I, respectively 
wave number q and the energy E is shown for the two cases of (a) weak disorder and (b) 
strong disorder. For both cases the sketch indicates the rough position of the different 
transport regimes, namely, localized (loc, gray), diffusive (dif., white), and ballistic 
(bal., black). 

It is well known that in the limit of weak disorder the localization phenomenon is 
restricted to the borders of the spectrum solely. Deep in the outer tails of the spectrum 
the states are localized on a single lattice site, whereas the overwhelming majority of 
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all states, i.e., not only the ones from the spectral middle, is still extended. Thus, as 
displayed in figure [U the localization length (envelope of gray areas) is not a closed curve 
in the (q, _E)-space. The energies which separate localized and non-localized regimes at 
q = (points between gray and white areas at q — 0) are the mobility edges. And in 
fact, much work has been devoted to the concrete position of the mobility edges [21 [3]. 
Only for energies between the mobility edges there is a conductor at the infinite length 
scale, otherwise there is an insulator at that length scale, of course. But, as already 
mentioned before, insulating behavior is practically absent in the limit of weak disorder, 
e.g., the Fermi level is much larger than the lower mobility edge. 

While insulating behavior appears at rather large length scales above the localization 
length, ballistic behavior occurs at comparatively small length scales below the mean 
free path (envelope of black area). Here, (quasi-)particles are not scattered and move 
freely with mean (group) velocities, e.g., as routinely evaluated in the framework of 
standard solid state theory. The latter free motion is reflected in the term ballistic and 
is typical for an ideal conductor. The mean free path, as drawn in figured], appears to be 
a contra-intuitive curve in the (q, £')-space, since it is smaller for states from the spectral 
middle than for states from the borders of the spectrum. However, we demonstrate in 
section 13.11 that in the limit of weak disorder such a curve results from the Boltzmann 
equation. 

Apparently, for weak disorder there is the practically unbounded regime above the mean 
free path where transport is neither insulating nor ballistic. This is the regime where 
transport is generally expected to be diffusive, i.e., at that length scales one expects a 
normal conductor. Particularly, figure [1] marks a "corridor" of wave numbers (dashed 
lines) with diffusive dynamics at almost all energies. (The notion of a diffusive corridor 
becomes helpful for later argumentations in the context of strong disorder where the 
existence of such a corridor is anything else than obvious.) From a mere theoretical 
point of view it is a challenge to concisely show that the dynamics is in full accord with 
a diffusion equation. But in section 13.11 we demonstrate by the use of the Boltzmann 
equation that the dynamics is indeed diffusive and further evaluate quantitatively the 
diffusion coefficient as a function of energy. As indicated in figure [U its dependence on 
energy seems to be as contra-intuitive as the one of the mean free path and strongly 
differs in the details from the approximations according to a free electron gas. 
For strong disorder the localization phenomenon is much more pronounced. When the 
amount of disorder is increased, the localized regimes gradually expand towards small 
length scales and towards energies in the middle of the spectrum as well, see figure [TJ On 
the one hand the already non-extended states become localized on smaller and smaller 
length scales. On the other hand more and more of the before extended states become 
localized at all. Hence, the mobility edges move closer to each other and eventually 
meet, once the critical disorder is reached. Then all states are localized and transport 
at the infinite length scale vanishes completely, i.e., at all energies. Therefore much 
work has addressed the concrete evaluation of the critical disorder [21 El S] - However, 
even above the critical disorder, transport takes place below the localization length, of 
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course. 

It is a priori not clear whether or not the dynamics below the localization length is still 
in good agreement with a diffusion equation in the limit of strong disorder, both below 
and above the critical disorder. But by the use of a method which is based on the TCL 
projection operator technique we demonstrate in section 13.21 that there also exists a 
corridor of wave numbers where the dynamics is indeed diffusive at almost all energies, 
at least as long as the amount of disorder does not become too strong. In particular 
the diffusion constant within this corridor does not substantially depend on energy, 
see figure HJ In fact, only if the dynamics for a certain wave number is not governed 
by a significant energy dependence, the method makes a definite conclusion, otherwise 
no information results except for the strong energy dependence of the dynamics, e.g., 
the method can not distinguish between highly energy dependent diffusion coefficients 
and non-negligible localized contributions. However, once a diffusive corridor of wave 
numbers with a single diffusion constant is reliably detected, it is natural to assume 
that the diffusion coefficient does not change, when this corridor is left. (Per definition 
diffusion coefficients should not depend on the wave number). Or, in other words, we 
suggest that in the limit of strong disorder the dynamics in the whole diffusive regime 
is well described by an energy independent diffusion constant. This diffusion constant 
is quantitatively evaluated in section 13.21 as a function of disorder. 
For the case of strong disorder the TCL-based method additionally allows to characterize 
the ballistic regime, i.e., by the use of the method the mean free path and the mean 
velocity can be also evaluated. Similarly, these quantities are found to be approximately 
independent from energy. This observation suggests that in the limit of strong disorder 
the whole dynamics below the localization length is not governed by significant energy 
dependencies. Of course, in a sense this suggestion disagrees with the observations for 
the case of weak disorder. Nevertheless, in the sections 13.11 and 13.21 the disagreement is 
subsequently resolved, both qualitatively and quantitatively. 

3. Quantitative results 

3.1. Weak disorder: Boltzmann equation 

In the present section we are going to investigate the dynamics in the limit of weak 
disorder. In that limit there certainly is a large variety of different approaches which all 
treat the disorder as a small perturbation to the clean Hamiltonian. Here, we briefly 
review on one class of these approaches, namely, the mapping of the quantum dynamics 
onto Boltzmann equations [T3l [HI [151 EH]. Different approaches to such a map rely on 
different assumptions and/or approximation schemes which are not entirely free of their 
own subtleties. However, the particle velocities that eventually enter the Boltzmann 
equation are routinely taken from the clean (unperturbed) Hamiltonian. To this end 
the clean Hamiltonian has to be diagonalized at first. Routinely, this diagonalization 
can be done by the application of the Fourier transform. Then the Hamiltonian takes 
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on the form 



H = ^E 9 aU 9 , (2) 



where at, a- are creation, annihilation operators for (quasi-)particles with the wave 
vector q, i.e., = 2-7T ki/N, kj = 0, ... ,N — 1. The corresponding dispersion relation 



reads 



£ 9 -=£2(l-cosg 4 )^ |g| 2 . (3) 



=i 



(Now and in the following the indicated approximations hold true for sufficiently small 
| q\, respectively low energies and are well known from the free electron gas.) As long 
as disorder is absent, the (quasi-)particles are not scattered and may be said to move 
freely with the (group) velocities which are determined by the derivative of the dispersion 
relation, namely, 



VfEf \=2 X Y, sin2 & ~ 2 I <fl ■ (4) 
\ i=i 

Whithin such a Boltzmann equation framework disorder takes the role of a set of 
impurities from which the (quasi-)particles are scattered after, say, a mean free time 
t(E), respectively mean free path 1(E). Therefore disorder essentially gives raise to a 
linear collision term, i.e., a rate matrix which describes the transitions between different 
(quasi-) momentum eigenstates. Generally, diffusion coefficients may be computed based 
on the inverse of this rate matrix. However, since the disorder of the Anderson model 
(statistically) features full spherical symmetry, a relaxation time approximation turns 
out to be exact, even though the dispersion relation does not feature full spherical 
symmetry [16]. Following this approach, the diffusion coefficient may be cast into the 
basic form 

V(E) = U(E)l(E), l(E) = v(E)r(E), r(E) = — -L— , (5) 

where v(E) denotes a mean velocity which is obtained from an average over all q 
featuring a certain energy E, i.e., 



v(E) = (v f ) { fl E9=jS} K2VE. (6) 
Furthermore, p(E) expresses the density of states normalized to the volume, i.e., 



1 dZ(E) VE 

P{E) = JP ^E~ W 4^ (7) 
w.r.t. the clean Hamiltonian. As a first observation, diffusion coefficients and mean 
free paths are inversely proportional to the amount of disorder, at least within the 
Boltzmann equation approach at hand. Due to the above mentioned subtelties of the 
mapping itself it is hard to give a detailed estimate for the regime of its applicability. 
However, disorder should generally be substantially smaller than regular hopping, i.e, 
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a < 1. 

Inserting the approximations for low energies in 



© into dSj) yields 



V(E) 



8ir^E 



4:71 

1(E) » - . 



(7" 



(8) 



This result coincides with the one in [2] which is found therein by the use of Green's 
functions. However, in order to obtain the full energy dependencies of T>(E) and 1(E) 
we numerically evaluate (E} and ([7]) in figure |2] (a) and (b). Note that the evaluation 
can be done for very large lattices, e.g., N = 1000, since exact diagonalization is not 
involved. Obviously, the approximations for and ([7]) are valid for very low energies 
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Figure 2. Energy-dependencies in the limit of weak disorder: (a) density of states, 
(b) mean velocity, (c) mean free path as well as (d) diffusion constant (according 
to the Boltzmann equation). The numerical data (circles) is extracted from a cube 
with 1000 3 lattice sites. The known approximations for low energies (solid lines) are 
indicated for comparison. 



solely, i.e., in the outer tails of the density of states. The actual curves differ strongly in 
the details. As a consequence the curves for T>(E) and 1(E), as displayed in figure 0(c) 
and (d), show interesting features, too. Particularly, the maximum diffusion coefficient 
is not located at the middle of the spectrum (E = 0). Instead two distinct maxima are 
observed at positions which are closer to the borders of the spectrum (E 4.5). 
The curve for T>(E) seems to already indicate that the overall dynamics of all energy 
regimes, i.e., the dynamics at high temperatures can not be described as diffusive with 
a single diffusion coefficient. However, for a definite conclusion the T>(E)-cmve has to 
be weighted with the density of states p(E), of course. Therefore in figure [3] the relative 
number of states r is shown which contribute to a certain diffusion constant. In a sense 
r is again a density of states but now in the space of diffusion coefficients. Apparently, 
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Figure 3. The relative number of states r corresponding to a diffusion constant V 
(according to the Boltzmann equation), cf. Fig.[2l The bars of the histogram visualize 
contributions from energies \E\ > 3 (light-colored area) and from energies \E\ < 3 
(dark-colored area). 



the majority of all states corresponds to diffusion constants which are rather close to 
V « 2.5/a 2 , i.e., the value for E = 0. These states are also located around the middle 
of the spectrum, as indicated in figure [3j But there is a relevant number of states from 
the outer parts of the spectrum which contribute to larger values of T>. Remarkably, in 
these parts the number of states with smaller values of T> is negligible. However, figure [3] 
clearly demonstrates that the overall dynamics of all energy regimes is not diffusive with 
a single diffusion coefficient. 

The situation may change, when the limit of weak disorder is slightly left, i.e., when 
the above predictions of the Boltzmann equation begin to break down. Obviously, the 
breakdown begins for the states from the borders of the spectrum, since these states 
are the first which become eventually localized. (The Boltzmann equation does simply 
not predict localization). At this point the predictions for the states from the spectral 
middle are still unaffected, of course. On that account it may happen that the large 
values of T> in figure |3] are gradually moved towards T> w 2.5/ a 2 such that r finally 
becomes more or less peaked at this position. In that case the dynamics is governed by 
a single diffusion constant. So far, this line of thoughts is a mere assumption. Even if 
the assumption was correct, it would be entirely unclear whether or not this assumption 
has some impact on a situation with strong disorder, i.e., beyond any validity of the 
Boltzmann equation. 

In the following sections 13.21 and 13.31 we subsequently show for strong disorder that it 
appears to be indeed justified to describe the dynamics below the localization length as 
diffusive with a single diffusion coefficient. Surprisingly, this diffusion constant is rather 
close to the value 2.5/a 2 , wide outside the strict validity of the Boltzmann equation. 
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3.2. Strong disorder: TCL projection operator technique 



Our approach in the limit of strong disorder is based on the time-convolutionless (TCL) 
projection operator technique [TFl [18] which has already been applied to the transport 
properties of similar models but without disorder, see [19j [20j [21] . I n its standard form 
this approach is restricted to the infinite temperature limit. This limitation implies that 
energy dependencies are not resolved, i.e., our results are to be interpreted as results on 
an overall behavior of all energy regimes. 

As illustrated in figure HJ we consider a 3-dimensional lattice consisting of N layers 
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Figure 4. A 3-dimensional lattice which consists of N layers with nxn sites each. Only 
next-neighbor hoppings are taken into account. Constants for intra-layer hoppings are 
set to 1 (white arrows), inter-layer hoppings are specified by another constant A (black 
arrows). 



with nxn sites each. The Hamiltonian of our model is almost identical to ([I]) with one 
single exception: All hopping terms which correspond to hoppings between layers (black 
arrows in figure HJ) are multiplied by some constant A. This multiplication is basically 
done due to technical reasons, see below. However, for A = 1 the Hamiltonian reduces 
to the usual Anderson Hamiltonian ([1]). 

We now establish a coarse-grained description in terms of sub units (similarly to |22j): 
At first we take all those terms of the Hamiltonian which only contain the sites of the 
fith layer in order to form the local Hamiltonian of the subunit //. Thereafter all 
those terms which contain the sites of adjacent layers \i and fi + 1 are taken in order 
to form the interaction Xv^ between neighboring subunits fi and fi + 1. Then the total 
Hamiltonian may be also written as 

N-l N-l 

H = H + \V, # =EV V=Y,*n ( 9 ) 

At=0 fi=0 

where we use periodic boundary conditions, e.g., we identify ji = N with \i = 0. The 
introduction of the additional parameter A allows for the independent adjustment of the 
interaction strength in this coarse-grained description. Since we are going to work in the 
Dirac picture, the indispensable eigenbasis of Hq may be found from the diagonalization 
of disconnected layers. 

By we denote the particle number operator of the /xth subunit, i.e., the sum of a\a ? 
over all f of the /xth layer. Because the overall number of particles is conserved, i.e., 
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[Y^nPut H\ = and no particle-particle interactions are taken into account, we still 
restrict the investigation to the one-particle subspace. The actual state of the system 
is naturally represented by a time-dependent density matrix pit), i.e., the quantity 
Pfj,{t) = Tr{ p(t) p M } is the probability for locating the particle somewhere within the 
pth subunit. The consideration of these coarse-grained probabilities corresponds to the 
analysis of transport along the direction which is perpendicular to the layers. Instead 
of simply characterizing whether or not there is transport at all, we analyze the full 
dynamics of the p^if). 

The dynamical behavior of the p^it) may be called diffusive, if the p^it) fulfill a discrete 
diffusion equation 

p,{t) = V \p^(t) - 2p M (t) + p, +1 {t) ] (10) 

with some p- and i-independent diffusion constant T>. It is a straightforward manner to 
show (multiplying (flUj) by p, respectively p 2 , performing a sum over p and manipulating 
indices on the r.h.s.) that the spatial variance 

2 

(11) 



N-l 

Var(t) = 

fi=0 



N-l 

Writ) 



increases linearly with t, i.e., Var(t) = 2T>t. Contrary, ballistic behavior is characterized 
by Var(i) oc t 2 , while insulating behavior corresponds to Var(i) = const., of course. 
According to Fourier's work, diffusions equations are routinely decoupled with respect 
to, e.g., cosine-shaped spatial density profiles 

P q (t) = c q J2 C0S (<?^)Pm(X) , ^ = ~7T' fc = M,---,y ( 12 ) 

and a yet arbitrary normalization constant C q . Consequently, (flOj) yields 

p q {t) = -2 {1- cos q)V Pq {t). (13) 

Therefore, if the quantum model indeed shows diffusive transport, all modes p q (t) have 
to relax exponentially. If, however, the modes p q {t) are found to relax exponentially only 
for some regime of q, the model is said to behave diffusively on the corresponding length 
scale / = 7c/q. One might think of a length scale which is both large compared to the 
mean free path (below that ballistic behavior occurs, <J 2 {t) oc t 2 ) and small compared 
the localization length (beyond that insulating behavior appears, <J 2 (t) = const.). 
For our purposes, i.e., for the application of the TCL projection operator technique, it 
is convenient to express the modes p q {t) as the expectation values of respective mode 
operators p q , namely, 

N-l 

p q (t) = Tr{ p(t)p q } , p q = C q cos(q p) , (14) 

fi=0 

where the normalization constants C q are now chosen such that Tr{p 2 } — 1. With this 
normalization 

V p{t) = Tr{ p{t) p q }p q = p q {t) p q (15) 
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defines a suitable projection (super)operator, because V 2 = V . For those initial states 
p(0) which satisfy V p(0) = p(0), i.e., for harmonic density profiles the TCL projection 
operator technique eventually leads to a differential equation of the form 

p q (t) = R q (t) Pq (t) , R q (t) = X 2 R 2 ,q{t) + X 4 R 4 , q (t) + . . . (16) 

which is a formally exact description for the dynamics at high temperatures, since p(0) 
is not restricted to any energy subspaces. Apparently, the dynamics of p q (t) is controlled 
by a time-dependent decay rate R q (t). This decay rate is given in terms of a systematic 
perturbation expansion in powers of the inter-layer coupling. (Concretely, for this model 
all odd orders vanish.) At first we concentrate on the truncation of (fTBT) to lowest order, 
i.e., to second order. But the fourth order is considered afterwards in order to estimate 
the validity of this second order truncation. 

According to [18], the TCL formalism routinely yields the second order prediction 

p q (t) = X 2 R 2 , q (t) Pq (t) , R 2>q (t) = /*dA/,(A) (17) 

Jo 

with the two-point correlation function 

f q (A) = Ti{[V(t),p q ][V(t'),p q }}, A = t-t', (18) 

where the time dependencies of operators are to be understood with respect to the 
Dirac picture. The g-dependence in (TT8T) is significantly simplified under the following 
assumption: The autocorrelation functions Tr{ v^t) v^it') } of the local interactions £> M 
should depend only negligibly on the layer number p (at relevant time scales). In fact, 
numerics indicate that this assumption is well fulfilled (for the values of a which are 
discussed here), once the layer sizes exceed ca. 30 x 30. Therefore first investigations 
may be based on the consideration of an arbitrarily chosen junction of two layers. The 
local interaction between these representative layers may be called Vq. The use of the 
above assumption simplifies ( [TBI to 

/ ff (A)w-2(l-cosg)/(A), /(A) = Tr{ v (t) v (t') } , (19) 

n 

where the g-dependence enters solely as an overall scaling factor [21 J. As a consequence 
the second order prediction at high temperatures reads 

p q {t) ^ -2(1- cos q)X 2 R 2 (t) Pq (t), R 2 {t)= /*dA/(A). (20) 

Jo 

This equation is already very similar to (TT3|) but still contains a time-dependent diffusion 
coefficient T>(t) = X 2 R 2 (t). However, it numerically turns out that /(A) behaves like a 
standard correlation function, i.e., it decays completely within some time scale tc- After 
this correlation time /(A) approximately remains zero and R 2 (t) takes on a constant 
value R 2 , the area under the initial peak of /(A). Numerics indicates that neither Te- 
nor R 2 depend substantially on n (at least for n > 30) such that both tc and R 2 are 
essentially functions of o. Since the correlation time tq apparently is independent from 
q and A, it is always possible to realize a relaxation time 



2(1 - cos q)X 2 R 2 
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which is much larger than tc, e.g., in an infinitely large system there definitely is a 
small enough q. For tr ^> Tc the second order prediction fl20|) at high temperatures 
immediately becomes 

p q (t) ^-2(1- cos q)X 2 R 2 p q (t) (22) 

and the comparison with ({TBI clearly shows diffusive behavior with a diffusion constant 
T> = A 2 i?2- Due to the independence of R2 from n (again for n > 30) the pertinent 
diffusion constant for arbitrarily large systems may be quantitatively inferred from a 
finite, e.g., 30 x 30 layer, see figure El Therein T> is evaluated for the range of a where 
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Figure 5. Theoretical prediction of lowest order TCL for an energy-independent 
diffusion constant T> as a function of the disorder a (Gaussian distribution) for the 
layer sizes n = 30 (crosses) and n = 90 (circles). The theoretical prediction is shown 
up to the critical disorder oc ~ 6 and the validity range of this prediction is displayed 
(dashed lines). As a guide for the eyes a proportionality to er~ 2 is indicated (solid 
line). For comparison the value T> = 1.05 ± 0.10 from [9] is shown (triangle), as found 
therein for E = and a = 2 (Gaussian distribution). Further values for V from [8] 
are displayed (triangles), as obtained therein for E — but for W = 1, 2 and 4 (box 
distribution) . In that case data for Gaussian and box distributions are supposed to be 
simply convertible by W/a w 2.6 (according to the ratio of the critical values). The 
prediction for E — according to the Boltzmann equation is also indicated (dotted 
line) . 



the used approximation for the g-dependence of the correlation function turns out to be 
justified, cf. f ll~8|) and ffl~9~]) . We additionally indicate already the validity range of the 
second order prediction at high temperatures, although this point is firstly discussed 
in detail in the next section 13.31 However, within the validity range there indeed is a 
corridor of q in which the dynamics at high temperatures can be described as diffusive 
in terms of (|22l) . Outside the validity range such a g-corridor does not exist, since either 
diffusion constants become highly energy dependent (a < 0.2) or localized contributions 
become non-negligible (a > 2), cf. figured! 

As indicated in figure |5l at the l.h.s. of the validity range the diffusion coefficient simply 
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scales as T> oc 1/cr 2 . Because such a scaling is expected in the limit of weak disorder, 
we quantitatively compare with the result 2.5/cx 2 for E — (A = 1) according to the 
Boltzmann equation, although the weak disorder limit is reached by no means, even 
for the energy regime around E = 0. The surprisingly good agreement supports the 
line of thoughts in section 13.11 namely, the diffusion constants at all energies are rather 
close to the value 2.5/a 2 , when the disorder becomes non-weak. At the r.h.s. of the 
validity range the diffusion coefficient begins to deviate from a simple l/<T 2 -scaling, e.g., 
D as 1.05 A 2 for a = 2. This value excellently agrees with the one in [9] which is found 
therein by a numerical study of Green's functions. The study in [9] remarkably requires 
an ensemble average over very many realizations of disorder, while a single disorder 
realization is adequate here, i.e., the correlation function is a self-averaging object. 
So far, we have characterized the diffusive regime. We now turn towards an investigation 
of the ballistic regime. Obviously, the replacement of ( 1201) by f l22l) is only self-consistent 
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Figure 6. Theoretical prediction of lowest order TCL for an energy-independent (a) 
mean free time r and (b) mean free path / w.r.t. the disorder a (Gaussian distribution) 
for the layer sizes n — 30 (crosses) as well as n = 90 (circles). As a guide for the eyes 
proportionalities to cr~ 2 are shown (solid lines). The predictions for E = according 
to the Boltzmann equation are also shown (dotted lines). The results in (a) and (b) 
suggest a mean velocity v/\ rj 2.5. 
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for tr ^> To, i.e., if the relaxation time is much larger than the correlation time. But 
this criterion has to break down for sufficiently large q, see figured! Hence, a transition 
from the diffusive to the ballistic regime appears for those modes p q (t) which decay 
on an intermediate time scale tr ~ tq, while the ballistic regime is finally reached for 
those modes p g (t) which decay on a short time scale tr <C tq. In fact, the rate i? 2 (£) 
is found to increase linearly for t < Tc/2ir. Due to this linear increase the second order 
prediction at high temperatures is a Gaussian decay of the corresponding modes. Such 
a Gaussian decay is known to be typical for ballistic dynamics [T9l [20], [2Tj . However, the 
ballistic character of the dynamics is most convincingly demonstrated in terms of the 
variance Var(t), as defined in (TTTT) . The second order prediction at high temperatures 
for this variance is given by 

Var 2 (t) = 2 A 2 f ' dt' R 2 (t') (23) 
Jo 

and scales as Var 2 (t) oc t 2 for t < t c /2tt. Of course, this result suggests the mean free 
time r = tc/2tt. Consequently, in order to obtain the mean free path as well, a mode 
p q (t) with tr = t has to be considered, i.e., the condition p q (r) /p q (0) = 1/e has to be 
fulfilled. This condition and (f20l) yield 

2(1 -cosg)A 2 [ T dt'R 2 (t') = 1 (24) 
Jo 

which can be rewritten as (1 — cosg) Var 2 (r) = 1. The use of I = ir/q eventually leads 
to the expression for the mean free path 

I 7T 

arccosfl — 1/Var 2 (r)] 

or the approximation / f» 71 / 'y/2 y / Var 2 (r) for I ^> 1, i.e., about two times the standard 
deviation. We use this approximation, because the dependence on A becomes trivial, 
namely, I oc A. However, whenever isslor even smaller, the concrete value of the mean 
free path is less important, since then the ballistic regime is restricted to a length scale 
below a single lattice site and does simply not exist. 

As indicated in figure El both the mean free time r and the mean free path / are 
proportional to 1/cr 2 over the full range of accessible a where the used approximation 
for the g-dependence of the correlation function turns out to be justified, cf. (fT8l) and 
( [191) . Therefore the mean velocity v — l/r becomes independent from a. Again there 
is a quantitative agreement with the prediction for E = (A = 1) according to the 
Boltzmann equation, wide outside the weak disorder limit. In contrast to figure El the 
validity range of the second order prediction at high temperatures is not indicated in 
figure E], since this prediction for short times is expected to be valid for all accessible 
a, see the next section [331 However, for o > 1.5 the mean free path / takes on values 
which are smaller than A, e.g., for A = 1 the ballistic regime is practically absent. 

3.3. Validity range of the TCL-based theory 

In the present section we are going to discuss the validity range of the second order 
prediction at high temperatures in more detail. To this end we consider the ratio 7Z(t) 
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of the fourth order to the second order, namely, 

R W = 2A-(S)ftW (26) 

cf. (|T6|) . Whenever 1Z(t) <C 1, the second order term 2A 2 (1 — cosg) R2, q (t) dominates 
the decay of the modes p q (t) and the fourth order term \ 4 R^ q (t) is negligible. But in 
general already the direct evaluation of the fourth order term turns out be extremely 
difficult, both analytically and numerically However, by the use of the techniques in 
the fourth order term can be approximated by 



A 4 R± :q {t) ^ 4 A 4 (1 - cos g) 2 i? 4 (t) (27) 
with the remaining g-independent rate 



R A {t) = t 



n 



\ E (/J dA (V>t| V*) Mt) l^)) 2 " R2(tf 



(28) 



where \ipi) denote the eigenstates of H . Consequently, in complete analogy to the rate 
R2{t), also the rate R^t) may be evaluated from the consideration of an arbitrarily 
chosen junction of two layers. The local interaction between these representative layers 
is still called £>o. The above approximation is based on the fact that the interaction V 
features the so-called Van Hove structure [25], [26], i.e., V 2 essentially is a diagonal matrix 
(in the eigenbasis of Hq). However, for the concrete derivation of this approximation we 
refer to [21] and concentrate on the implications here. By the use of the approximation 
the ratio 7Z(t) can be rewritten as 

K(t)^Q^, Q = 2A 2 (l-cosg). (29) 

This ratio is a monotonically increasing function of t, cf. ( )28|) . As a consequence there 
always exists a time t# with lZ(tB) = 1, i.e., a time where the contributions R2(t) and 
i?4(t) are equally large. But this fact does not restrict the validity of the second order 
prediction, if t# ^> tr and hence 1Z{tr) 1. The validity obviously breaks down only 
in the case of, say, IZijn) ps 1 or even larger. Since both 7Z(t) and tr depend on Q, we 
use again the condition p q {r R ) / p q {Q) = 1/e, i.e., 



Q dt'R 2 {t') = l (30) 
Jo 

in order to replace Q in fl29l) . Due to this replacement 1Z(tr) becomes a function 

K(Te) = R 2 (r R) R $i R 2 (t) (31) 

of the free variable tr. (tr still depends on Q, of course.) Because also 1Z(tr) increases 
monotonically, we define max(r^) as the maximum tr for which 1Z(tr) is still smaller 
than 1. This maximum relaxation time already specifies the validity range of the second 
order prediction. However, it is useful to set max(r^) in relation to the correlation time 
t c . We therefore define the measure x as the dimensionless quantity x — r c/ max ( r ij)) 
e.g., x = 1 directly implies the breakdown of the second order prediction on relatively 
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Figure 7. The measure \ f° r the validity range of lowest order TCL as a function of 
the disorder a and the inverse layer size 1/n 2 . This measure detects the "corridor" of 
wave numbers, where transport is diffusive at almost all energies with a single diffusion 
constant, cf. Fig. [TJ For those values of x which are on the order of 1/4 the corridor 
does not exist. But for those values of \ which are closer to zero the corridor opens. 
The smaller x, the larger is this corridor. An absolute minimum Xmin ~ 1/50 is found 
at a m 1/2 in the limit n — > oo. Note that only 10% of the whole area is extrapolated 
(the area in front of the thick line) . 



short time scales on the order of tc, whereas x = strongly indicates its unrestricted 
validity. For practical purposes an interpretation of x hi the context of length scales 
certainly is advantageous. Such an interpretation essentially requires the inversion of 
f !3Up . In general this inversion can only be done by numerics. But we have tr = 1 /(Q.R2) 
for tr ^> Tc and may hence write 

a r~ = c ' o ■ p = — 9 — ' ( ' 

amax -K2 ^mm n 2 A 

where the factors 2 and 1/2 are chosen to slightly fulfill r c < r fl < max(rij), i.e., 
Qmax, Q m in correspond to lower, respectively upper border of the diffusive corridor, 
cf. figure [TJ We finally end up with 



^min 



4X (33) 



Qmax 

or by the use of Q ~ q 2 X 2 with g m j n /gmax ~ 2y/X- Thus, for those values of x which 
are on the order of 1/4 the corridor does not exist. But for those values of x which are 
closer to zero the corridor opens. The smaller \i the larger is this corridor. 
In figure [7] the measure x is quantitatively evaluated as a function of the amount of 
disorder a and the inverse layer size 1/n 2 . (The rate Ri(t), other than the rate R2{t), 
scales significantly with n. This scaling gives rise to the n- dependence of x-) F° r each 
layer size there is a optimum disorder where x is minimized, i.e., where the diffusive 
corridor is maximized. But for n = 30 (back of figure [7]) we find l^fx ~ 2/3 at the 
optimum disorder. This value indicates a corridor of about one or two diffusive modes 
(for N = 30). For all a and n < 100 (which is the limit for our numerics) x clearly 
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appears to be of the form 

X (a,n) = ^ + B(a). (34) 

The extrapolation of the 1 / n 2 -scaling eventually leads to a suggestion for n = oo (front 
of figure [TJ. According to this suggestion, we find 2^/x ~ 2/7, again at the optimum 
disorder. This value indicates a still narrow but existent corridor of diffusive modes (for 
N = oo). 

We finally recall that these findings apply at infinite temperature, i.e., the narrow 
diffusive corridor is characterized by the fact that the dynamics within this corridor 
is diffusive at almost all energies with a single diffusion coefficient. The narrowness of 
this corridor passes into a complete absence, since either diffusion constants become 
highly energy dependent (er < 0.2) or localized contributions become non-negligible 
(a > 2), cf. figure [H 



3.4- Numerical verification 

In the last two sections we have introduced the TCL-based method and have discussed 
its predictions as well as the validity of these predictions. In the present section we are 
going to present the results of numerical simulations in order to verify the predictions 
of the method, as far as possible from the consideration of a finite system. Since the 
applicability of the method requires a system which consists of layers with a minimum 
size of n = 30, we consider layers of that size in the following simulations. According 
to the predictions of the method, for n = 30 a diffusive corridor is only existent for 
disorders in the vicinity of a — 1, see figure UJ Therefore we focus on such a value of a 
in all numerical simulations. 

For cx=l the TCL-based theory predicts a diffusion constant T> « 2.9 A 2 and a mean 
free time r ~ 1.1, i.e., a correlation time tq = 2ht 6.9, cf. figures [5] and [6j According 
to the theory, diffusive dynamics emerges only on a time scale which is given by the 
condition tq <C Tr = l/[2 (1 — cosg) T>], e.g., a sufficiently small q has to be chosen. For 
the naturally interesting isotropic case of A = 1 the choice q = n/20 leads to the ratio 
tr/tc ~ 2.0. Unfortunately, q = n/20 is firstly realized for a system which consists of 
N = 40 layers and such a system already is too large for the application of numerically 
exact diagonalization. 

However, approximative numerical integrators may be applied, e.g., on the basis of a 
Suzuki- Trotter decomposition of the time evolution operator [271 1211 [29]. In detail we 
choose a pure initial state \ijj q (0)) and apply a fourth order Suzuki- Trotter integrator in 
order to obtain the time evolution \ip q {t)) of this initial state and to evaluate the actual 
expectation value p q (t) = {%fj q {t)\p q \%fj q {t)). In particular we choose the initial state at 
random and only require the condition (ip q (0)\p q i \ip q (0)) = S qiq ', i.e., we still consider a 
harmonic density profile. The result of the approximative numerical integrator is shown 
in figure El for a single realization of \ip q (Q)) with q = tt/20. Apparently, there is a very 
good agreement between this result and the prediction of the TCL-based theory. The 
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time t 

Figure 8. Results for the time evolution of a mode p q (t) with q = 7r/20 (without 
any restriction to energy regimes) . The numerical data (circles) is obtained by the use 
of a 4th order Suzuki- Trotter integrator for the model parameters N = 40, n = 30, 
(7 = 1 (Gaussian distribution) and A = 1 (isotropic hopping constants). The theoretical 
prediction of lowest order TCL (solid line) is indicated for comparison. 

latter agreement further demonstrates that the validity of the theoretical prediction is 
not restricted to an initial density matrix of the strict form p(0) = p q . This fact may 
be understood in terms of dynamical typicality [301 EH [32j [331 EI]- 
Although the above Suzuki- Trotter integrator allows to determine the time evolution 
of pure initial states for rather large systems, this integrator is not able to resolve the 
energy dependencies of the dynamics, of course. To this end we have to use numerically 
exact diagonalization which is applicable to a maximum system with about N = 10 
layers. In such a system diffusive dynamics is expected to emerge only, if the coupling 
constant A is much decreased. For A 1, in complete analogy to the above numerical 
simulation, the time evolution of pure initial states has comprehensively been shown 
to be in full accord with all predictions of the TCL-based theory [211 GO]- However, 
since it still remains to resolve the energy dependencies of the dynamics, we consider 
the quantities 

p q>E (t)=Tr{P E p(t)P E p q }, (35) 

where P E denotes a projector onto the states of some energy regime E. In pratice we 
choose a coarse-grained partition into five energy intervals with the same number of 
states, namely, there are 1800 states in each energy interval. A fine-grained partition 
into more energy intervals is not convenient, because only a sufficiently coarse-grained 
partition assures p g (t) ~ J2 E Pq,E{t), i.e, the quantities p q , E (t) resolve the dynamics on 
a reasonable energy scale. (The longer the relevant time scale for the dynamics, the 
smaller is this reasonable energy scale, of course.) 

The quantities p q>E (t) can be used in order to provide a diagram for the dependence 
of transport on q and E, similar to the sketch in figure [TJ For instance, two measures 
for the deviation of p q , E (t) from a strictly exponential decay (diffusive behavior) may 
be defined, cf. [35]: The first measure detects deviations towards a Gaussian decay 
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Figure 9. Numerical results for the deviation of transport from diffusive towards (I) 
localized and (II) ballistic behavior as a function of the wave number q and the energy 
interval i for inter-layer hopping constants (a) A = 0.08, (b) A = 0.09 and (c) A = 0.16. 
The color palette is chosen from white (small deviations) to black (large deviations). 
The underlying data is obtained from exact diagonalization for the model parameters 
N = 10, n = 30 and a = 1 (Gaussian distribution). 

at short times (ballistic behavior) and the second measure detects deviations towards 
a stagnant decay at long times (localized behavior). Such measures are displayed in 
figure El Whenever one of these measures is large (black areas), p q ,E(t) does not relax 
exponentially. But whenever both measures are small (white areas), p g ,E(t) decays 
exponentially, i.e., it behaves diffusively. Particularly, there indeed is a g-corridor where 
Pq,E{t) decays exponentially for practically all E. As predicted by the TCL-based theory, 
this diffusive corridor is shifted to smaller q, when A is increased. According to figure [9] 
(c), the borders q m [ n , qmax of the diffusive corridor lead to a ratio q m [ n / qm&x between 
1/2 and 1. The latter ratio remarkably is in accord with the theoretical prediction 
g m j n /gmax ~ 2/3, too. (In general the ratio Q m j n /<2max has to be compared. But 
for the borders in figure |9] (c) the approximation Q w A 2 q 2 is already justified.) 
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However, it still remains to clarify whether or not the dynamics within the diffusive 
corridor is governed by a single diffusion coefficient. To this end an exponential fit may 
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Figure 10. Results for the diffusion coefficient T> as a function of the inter-layer 
hopping constant A for the model parameters n — 30 and a = 1 (Gaussian distribution) . 
The data corresponding to figure [9] is indicated by crosses (mean value of energy 
intervals) with bars (mean deviation of energy intervals). Additional data without 
any restriction to energy regimes is presented for N — 10 from exact diagonalization 
(circle) and for N = 10, 20, 30, 40 by the use of a 4th order Suzuki- Trotter integrator 
(triangles), cf. figure M 

be applied to Pq^if) for each E in this g-corridor. Such a fit directly yields a decay 
rate and consequently a diffusion coefficient T>e- Then the mean value T> and the mean 
deviation 5T> of the energy intervals may be evaluated, see figure [10J For all A where a 
diffusive corridor exists for iV = 10 we find a mean value T> between 2.9 A 2 and 3.0 A 2 
and a mean deviation 5T> on the order of less than 10%. This finding finally supports 
the theoretical prediction for a diffusive corridor with a single diffusion coefficient. For 
completeness, figure [10] additionally shows diffusion coefficients from a exponential fit 
to p q (t), as obtained by the use of the above Suzuki- Trotter integrator for pure initial 
states. Even though 5T> is not availabe in that case, T> scales simply as T> oc A 2 for all 
A up to the isotropic case of A = 1. 

4. Summary and conclusion 

In the work at hand we have investigated single-particle transport in the 3-dimensional 
Anderson model with Gaussian on-site disorder. Particularly, our investigation has been 
focused on the dynamics on scales below the localization length. The dynamics on those 
scales has been analyzed with respect to its dependence on the amount of disorder and 
the energy interval. This analysis has especially included the quantitative evaluation of 
the characteristic transport quantities, e.g., the mean free path which separates ballistic 
and diffusive transport regimes. For these regimes mean velocities, respectively diffusion 
coefficients have been evaluated quantitatively, too. 

By the use of the Boltzmann equation in the limit of weak disorder we have shown that 
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all transport quantities substantially depend on the energy interval. In addition we have 
demonstrated that these energy dependencies significantly differ from the well known 
approximations for a free electron gas. This significant difference develops for energies 
around the spectral middle where the overwhelming majority of all states is located. As 
a consequence the diffusion coefficients for these energies seem to be both a new and 
relevant result. 

In the limit of strong disorder we have found evidence for much less pronounced energy 
dependencies by an application of a method on the basis of the TCL projection operator 
technique. This method suggests that all transport quantities take on values which are 
practically independent from the energy interval. Remarkably, the latter values coincide 
with the prediction of the Boltzmann equation for the spectral middle, if this prediction 
is simply extrapolated to strong disorders, i.e., to disorders beyond any strict validity of 
the Boltzmann equation. Solely the suggested diffusion coefficient begins to differ from 
such a simple extrapolation, once the amount of disorder becomes on the order of the 
critical disorder. In the strict sense the TCL-based method does not yield a diffusion 
constant in the close vicinity of the critical disorder, because the validity range of the 
method is left for such an amount of disorder. In the close vicinity of the critical disorder 
the diffusion constant has to be understood as a mere conjecture. However, the method 
leads to a reliable diffusion coefficient for strong disorders which pass through almost 
one order of magnitude. Such a comprehensive description appears to be novel in the 
literature. 

Strictly speaking, the TCL-based theory makes only a definite conclusion on a corridor 
of finite length scales where the dynamics is diffusive at approximately all energies with 
a single diffusion coefficient. But we do not expect that the diffusion coefficient in the 
diffusive regime outside this corridor is significantly different, especially since diffusion 
constants should not depend on the length scale per definition. The latter expectation is 
also supported by the agreement with the Boltzmann equation and with the numerical 
results for diffusion constants in J8j |9]. However, whenever the above corridor of length 
scales is not existent, the theory does not allow for any conclusion. Since such a corridor 
may not exist in lower dimensions, the TCL-based theory may not lead to results on 
transport in the one- or two-dimensional Anderson model. But the theory itself, as 
demonstrated for the three-dimensional case, can analogously be applied also to the 
lower-dimensional cases, of course. This application is a scheduled project for the near 
future. 
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